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Predictions  are  reported  of  the  size  response  of  various  light-scattering 
aerosol  counters  manufactured  by  Particle  Measuring  Systems.  Models 
considered  are  those  that  exploit  the  high  intensity  of  light  available  within 
the  cavity  of  a  HeNe  gas  laser  (generically  referred  to  by  the  manufac¬ 
turer  as  "active  scattering  aerosol  spectrometer  probes").  The  new  response 
function  properly  averages  over  particle  trajectories  through  nodes,  anti¬ 
nodes,  and  intermediate  regions  of  the  intracavity  laser  beam.  Our  studies 
address  probes  having  two  basic  scattering  geometries:  those  that  collect 
light  scattered  over  a  relatively  narrow  solid  angle  (subtending  angles  be¬ 
tween  4°  and  22°  from  the  laser  beam  axis)  and  those  that  collect  light  over 
a  rather  large  solid  angle  (between  35°  and  120°).  The  new  response  func¬ 
tion  predicts  smoother  dependence  on  particle  size  than  the  previous  re¬ 
sponse  function  of  Pinnick  and  Auvermann  (1979,  J.  Aerosol  Sci.  10: 55-74) 
and  is  in  better  agreement  with  measurement.  Response  calculations  for 
common  atmospheric  aerosol  (water,  sulfuric  acid,  ammonium  sulfate,  and 
black  carbon)  reveal  the  considerable  sensitivity  of  the  response  to  particle 
dielectric  properties.  Comparison  of  response  calculations  with  the  manu¬ 
facturer's  calibration  reveals  conditions  for  which  the  manufacturer's  cali¬ 
bration  is  most  appropriate,  and  the  potential  for  errors  (as  much  as  a  factor 
of  two  in  sizing)  when  it  is  blindly  applied.  These  results  should  help  the 
user  of  these  instruments  to  more  realistically  interpret  size  distribution 
measurements. 
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1.  Introduction 


Light-scattering  particle  counters  are  among  the  most  widely  used  instru¬ 
ments  for  aerosol  measurement.  They  have  been  employed  for  determin¬ 
ing  estimates  of  the  tropospheric  and  stratospheric  aerosol  burden,  for 
monitoring  concentrations  of  particles  in  clean  rooms  essential  for  high- 
technology  manufacturing,  for  determining  aerosol  filter  efficiencies,  for 
detecting  atmospheric  aerosol  pollutants,  for  monitoring  particulates  in 
work  place  environments,  etc. 

These  devices  work  as  follows.  Air  containing  particles  is  drawn  through 
an  illuminated  volume,  where  light  scattered  by  single  particles  is  sensed 
and  converted  to  an  electrical  signal  whose  pulse  height  is  analyzed.  The 
pulse  height  is  used  to  infer  particle  size.  The  accumulated  measurement  of 
many  particles  results  in  a  size  distribution.  Particle  concentration  is  deter¬ 
mined  from  total  counts.  Optical  particle  counters  have  the  advantage  of 
providing  rapid,  nonintrusive  measurement  of  micrometer-sized  aerosol 
particles,  and  can  run  continuously  with  only  electrical  power  required  for 
their  operation. 

About  25  years  ago.  Particle  Measuring  Systems  (PMS)  in  Boulder,  CO,  be¬ 
gan  making  several  laser-based  light-scattering  particle  counters  (Schuster 
and  Knollenberg,  1972;  Knollenberg,  1976)  for  use  by  the  aerosol  commu¬ 
nity.  During  the  late  1970s  and  early  1980s,  a  number  of  these  instruments 
became  available.  The  manufacturer  refers  to  these  counters  as  laser  aerosol 
spectrometers,  with  the  implication  that  they  are  superior  to  other  light¬ 
scattering  particle  counters.  In  some  of  these  instruments  (the  ASASP-300, 
ASASP-300X,  ASASP-100X,  ASASP-X,  LAS-250X,  LAS-X,  HS-LAS,  a  num¬ 
ber  of  aircraft  versions  of  these  probes,  and  various  one-of-a-kind  models), 
particles  are  illuminated  by  the  intracavity  light  from  a  TEM0o  mode  HeNe 
laser.  The  manufacturer  generally  designates  these  probes  as  "active  cav¬ 
ity"  (probes  in  which  particles  are  illuminated  outside  the  cavity  are  des¬ 
ignated  "passive  cavity").  The  active-cavity  probes  are  the  subject  of  this 
report. 
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There  are  basically  two  scattering  geometries  used  in  the  active-cavity 
probes.  In  the  ASASP-300  and  ASASP-300X  probes,  the  scattered  light  is 
collected  by  a  lens  over  a  relatively  small  and  narrow  solid  angle  sub¬ 
tending  angles  between  4°  and  22°  from  the  laser  beam  axis.  In  the  more 
sensitive  ASASP-X,  ASASP-100X,  LAS-250X,  LAS-X,  and  HS-LAS  probes/ 
scattered  light  is  collected  by  a  parabolic  mirror,  allowing  collection  over 
a  relatively  large  and  wide  solid  angle  subtending  angles  between  35°  and 
120°  from  the  laser  beam  axis.  In  both  types  of  probes,  the  collection  solid 
angle  has  symmetry  about  the  laser  beam  axis.  The  "X"  in  the  model  desig¬ 
nation  denotes  that  the  probe  is  plumbed  with  a  relatively  small  inlet  orifice 
(typically  200  pm  in  diameter),  so  that  aill  particles  drawn  through  it  inter¬ 
sect  the  illuminating  laser.  Because  of  the  high  intensity  available  within 
the  cavity  (more  than  1  W  intracavity  circulating  power  for  most  models 
and  up  to  25  W  for  the  high-sensitivity  HS-LAS  model),  particles  as  small 
as  0.1  pm  in  diameter  can  be  measured. 

However,  even  with  their  novel  design  and  considerable  sensitivity,  these 
instruments  are  not  without  limitations:  in  particular,  their  size  resolution, 
which  we  address  in  this  report.  Their  sizing  limitations  arise  from  a  funda¬ 
mental  fact:  the  light  scattered  by  a  particle  depends  on  the  particle's  dielec¬ 
tric  properties  and  shape  as  well  as  its  size.  Thus,  unless  particles  of  known 
refractive  index  are  being  measured,  only  an  approximate  determination  of 
particle  size  can  be  made  from  the  magnitude  of  its  light-scattering  signal. 

Roughly  50  ASASP-300, 125  ASASP-X,  400  LAS-X,  and  60  HS-LAS  probes 
have  been  sold  by  PMS  to  date.  Because  of  their  wide  use,  the  performance 
of  these  instruments  has  been  studied  rather  extensively  (Pinnick  and  Au- 
verman,  1979;  Pinnick  and  Rosen,  1979;  Allan  and  Ashdown,  1982;  Garvey 
and  Pinnick,  1983;  Solderholm  and  Salzman,  1984;  Chen  et  al,  1984;  Liu  et 
al,  1985;  Szymanski  and  Liu,  1986;  Yamada  et  al,  1986;  Hinds  and  Kraske, 
1986;  Barnard  and  Harrison,  1988;  Knollenberg,  1989;  Jeung,  1990;  Kim  and 
Boatman,  1990;  Pueschel  et  al,  1990;  Hering  and  McMurry,  1991;  Liu  et  al, 
1992;  Knollenberg  and  Veal,  1992;  Kim,  1995). 

Many  of  these  studies  concentrate  on  sizing  performance  where  rather 
definitive  measurements  of  response  have  been  carried  out  on  carefully 
prepared  monodisperse  aerosol  test  particles  of  different  size,  shape,  and 
composition.  Other  studies  have  carried  out  theoretical  investigations  of 
particle  response  characteristics  using  the  approximate  response  function 

*Also  the  Royco  236  counter,  for  which  the  optical  system  was  manufactured  by  PMS 
but  sold  in  the  early  1980s  by  Hiac-Royco,  Menlo  Park,  CA. 


proposed  by  Pinnick  and  Auvermann  (1979)  and  Garvey  and  Pinnick 
(1983)  (eq  (2)  of  Pinnick  and  Auvermann,  1979,  hereinafter  referred  to  as 
the  PAG  response  function). 

For  the  narrow-angle  scattering  probes  (ASASP-300  and  ASASP-300X  mod¬ 
els),  these  studies  reveal  (1)  a  considerable  sensitivity  of  the  response  to 
particle  refractive  index;  (2)  the  existence  of  size  regions  where,  because  of 
Mie  scattering  resonances,  particles  of  different  size  give  the  same  response; 
and  (3)  degradation  in  size  resolution  for  nonspherical  particles  because 
of  the  dependence  of  scattering  on  particle  orientation.  The  same  findings 
generally  apply  to  the  more  commonly  used  wide-angle  scattering  probes 
(ASASP-X,  ASASP-100X,  LAS-X,  LAS-250X,  HS-LAS,  and  Royco  236  mod¬ 
els).  For  both  narrow-angle  and  wide-angle  scattering  probes,  measure¬ 
ments  of  the  response  in  the  cited  studies  generally  agree  with  the  PAG 
response  function,  although  the  high-frequency  oscillations  predicted  by 
the  PAG  response  function  were  not  found. 

Although  the  PAG  response  function  accounts  for  particle  index  of  refrac¬ 
tion,  it  is  only  approximate,  because  it  does  not  account  for  particles  travers¬ 
ing  nodes,  antinodes,  and  intermediate  regions  of  the  intracavity  standing 
wave  of  the  laser,  nor  does  it  account  for  laser  mode  changes  during  par¬ 
ticle  transit.  Averaging  over  particle  trajectories  and  phase  smooths  the  re¬ 
sponse  function  and  could  result  in  better  agreement  with  measurement. 

For  the  LAS-X  and  the  high-sensitivity  HS-LAS  probes,  there  may  be  ad¬ 
ditional  smoothing  caused  by  multiline  lasing.  In  these  probes,  according 
to  the  manufacturer  (Knollenberg,  1989),  the  laser  used  has  sufficient  gain 
to  be  above  threshold  for  lasing  in  multiple  lines.  The  PAG  response  func¬ 
tion  does  not  account  for  multiline  lasing.  Nor  does  it  treat  inhomogeneous 
particles. 

We  provide  here  an  improved  response  function,  which  correctly  aver¬ 
ages  over  particle  trajectory  through  the  cavity  of  a  laser  and  accounts  for 
the  laser  coherence  time  being  shorter  than  the  time  of  transit  of  particles 
through  the  laser  beam.  This  response  function  may  have  been  suggested 
previously  at  a  conference  (Soderholm  and  Salzman,  1984),  but  there  is  in¬ 
sufficient  information  in  the  conference  proceedings  to  judge.  The  new  re¬ 
sponse  function  was  suggested  previously  by  Garvey  and  Pinnick  (1983) 
(their  eq  (2))  but  not  recommended  by  them. 

We  compare  the  new  response  function,  which  is  a  smoother  function  of 
particle  size  than  the  PAG  function,  to  previous  response  measurements 
of  Pinnick  and  Auvermann  (1979),  Szymanski  and  Liu  (1986),  Garvey  and 
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Pinnick  (1983),  and  Chen  et  al  (1984).  The  new  response  function  is  in  bet¬ 
ter  agreement  with  measurements  of  test  aerosols  compared  to  the  PAG 
function.  Thus,  the  measurements  support  the  superiority  of  the  new  more 
general  response  function  over  the  PAG  response  function. 

Here  we  provide  response  calculations  for  both  narrow-angle  and  wide- 
angle  active-scattering  probes,  using  the  new  response  function,  for  both 
homogeneous  and  inhomogeneous  particles  characteristic  of  atmospheric 
and  pollutant  aerosol.  For  reference,  we  compare  the  results  to  the  manu¬ 
facturer's  calibration,  and  to  the  previous  PAG  response  function. 

Finally,  we  present  response  calculations  for  a  multiline  laser  source  that 
may  be  used  in  the  HS-LAS  or  LAS-X  probes.  Contrary  to  the  manufac¬ 
turer's  claim,  multiline  operation  has  little  effect  on  the  probe  response 
characteristics. 

The  results  should  enable  users  of  these  probes  to  better  assess  their  sizing 
limitations. 


2.  Theoretical  Response  Function 


Particles  detected  and  sized  in  the  active-cavity  scattering  probes  are  il¬ 
luminated  by  the  intracavity  light  from  a  TEMoo-mode  HeNe  laser.  The 
particle  cross  section  for  light  scattered  into  the  probe  collection  aperture 
(defined  as  the  probe  response  function  R)  should  account  for  two  averag¬ 
ing  effects.  First,  because  particles  do  not  pass  perpendicularly  through  the 
beam,  they  traverse  nodes,  antinodes,  and  intermediate  regions  (fringes) 
of  the  intracavity  standing  wave  of  the  laser.  Second,  the  coherence  time 
of  the  laser  is  generally  shorter  than  the  particle  transit  time,  causing  the 
fringe  positions  to  change  during  particle  transit.  In  appendix  A,  we  de¬ 
rive  an  expression  for  the  response  function  that  accounts  for  both  of  these 
effects.  For  spherical  homogeneous  particles,  the  result  reduces  to 


R 


7 r  /A'" 

J Or,,  i  n 


sin  9  d,9 


|Si(0)|2  +  |S2(0)| 2  +  \Si(n  -  9) |2  +  |S2(tt  -  6) |5 


(1) 


where  k  is  the  wavenumber,  6  is  the  scattering  angle,  and  S\{9)  and  S2(0) 
are  the  amplitude  scattering  matrixes  (Bohren  and  Pluffman,  1983). 

In  appendix  A  we  also  derive  a  corrected  expression  for  the  PAG  theory, 
which  assumes  that  (1)  the  particle  trajectory  is  perpendicular  to  the  in¬ 
tracavity  standing  wave,  (2)  the  particle  passes  through  an  antinode,  and 
(3)  the  coherence  time  of  the  laser  is  longer  than  the  particle  transit  time 
through  the  beam.  The  result  is 

/  sind0[|Si(0)  +  5i(tt  -  9) |2  +  |S2(0)  -  S2(vr  -  0)|2].  (2) 


This  response  function  is  slightly  different  from  the  PAG  response  function: 

qr  max  „  n 

sm0^[|51(0)  +  51(7r-0)|2  +  |52(0)  +  52(7r-^)|2].  (3) 

k  lOnnn 


Finally,  we  address  the  issue  of  multiline  lasing.  Lasers  in  some  of  the 
active-cavity  probes  have  high  gain  and  lase  on  multiple  lines.  We  show 
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in  appendix  A  that  for  multiwavelength  operation,  the  response  function 
is  a  sum  of  the  response  functions  at  each  wavelength  weighted  by  the  ap¬ 
propriate  irradiance  fraction  in  each  wavelength. 


3.  Response  Measurements 


A  definitive  test  of  the  validity  of  the  theoretical  response  function  (eq 
(1))  requires  measurement  of  well-characterized  monodisperse  test  parti¬ 
cles  having  a  range  of  sizes  and  refractive  indexes.  For  active-cavity  probes 
there  is  an  abundance  of  measurements  available.  Here  we  concentrate  on 
the  measurements  of  Pinnick  and  Auvermann  (1979)  for  the  ASASP-300, 
Garvey  and  Pinnick  (1983)  for  the  ASASP-X,  Chen  et  al  (1984)  for  the  Royco 
236,  and  Syzmanski  and  Liu  (1986)  for  the  ASASP-300X  and  LAS-X. 

We  first  consider  the  narrow-angle  scattering  probes  (ASASP-300  and 
ASASP-300X).  Figure  1  shows  the  measured  response  to  nonabsorbing 
polystyrene  latex  and  dioctylphthalate  (DOP)  particles  and  to  highly  ab¬ 
sorbing  nigrosin  dye  and  carbon  black  particles,  compared  to  the  theoreti¬ 
cal  response  function  (eq  (1)).  For  each  probe,  a  single  normalization  con¬ 
stant  relating  response  voltage  to  particle  cross  section  is  determined  from 
a  fit  of  the  experimental  data  for  polystyrene  latex  to  the  theoretical  curve 
for  latex.  This  fitting  procedure  is  described  in  appendix  B. 

The  measurements  confirm  the  predicted  oscillatory  character  of  the  re¬ 
sponse  curve  for  nonabsorbing  latex  and  DOP  particles,  although  the  fine 
structure  cannot  be  resolved.  The  measurements  also  generally  corroborate 
the  smooth  monotonic  response  predicted  for  highly  absorbing  particles. 
As  noted  previously  by  Pinnick  and  Auvermann  (1979),  the  fall-off  in  the 
observed  response  for  highly  absorbing  nigrosin  particles  for  radii  greater 
than  1  jum  is  believed  to  be  due  to  a  reduction  of  laser  power  caused  by 
the  particles.  The  disagreement  between  carbon  black  measurements  and 
theory  could  result  from  the  difficulty  in  preparing  homogeneous  spherical 
particles  (having  refractive  index  1.775  —  0.508i  as  used  in  the  theory)  from 
India  ink. 

A  similar  comparison  for  several  models  of  the  wide-angle  probes  (the 
ASASP-X,  LAS-X,  and  Royco  236)  is  shown  in  figure  2.  Again,  we  empha¬ 
size  that  we  determined  a  single  normalization  factor  for  each  probe  by 
doing  a  weighted  fit  of  the  polystyrene  latex  measurements  to  the  theoreti¬ 
cal  response.  The  theoretical  response  is  generally  within  the  error  of  mea¬ 
surement,  or  what  appear  to  be  errors  in  measurement,  as  evidenced  by 
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Figure  1.  Comparison 
of  theoretical  response 
(eq  (1))  and  measured 
response  for  PMS 
narrow-angle  scattering 
aerosol  probes.  Theory 
and  measurements  are 
shown  for 
nonabsorbing 
polystyrene  latex  and 
dioctylphthalate,  and 
highly  absorbing 
nigrosin  dye  and 
carbon  black.  Data  sets 
are  offset  by  1  to  3  ^ 

orders  of  magnitude  on  -§ 
vertical  scale,  as 

co 

indicated.  A  single  Q- 

normalization  factor  for  £ 
each  probe  relates  ^ 

experimental  response  w 
voltage  to  o 

particle-scattering  cross  <g 
section.  Predicted  broad  CC 
Mie  scattering 
resonance  features  are 
corroborated  by 
measurements. 


4°  to  22°  collection 


Theory 


PSL  /(1 .592-0 0 


DOP/(1 .4845  -  0/) 


iMigrosin  dye 
(1.67-0.26/1^. 

O 

O^carbon  black 
(1.775-0.508/) 


f  Measurements: 

OOOOO  Szymanski  &  Liu  (1986)/ASASP-300X 
1 .0  V  =  1 .0  x  1 0-8  cm2 

•  •  •  •  •  Pinnick  &  Auvermann  (1 979)/ASASP-300 
1.0  V  =  5.8  x  10“9  cm2 


5  6  7  8  9 

10"1 


5  6  7  8  9 


Particle  radius  (pm) 


Figure  2.  Same  as 
figure  1,  except  for  PMS 
wide-angle  scattering 
aerosol  probes. 
Predictions  agree  with 
measurements  and 
reveal  a  generally 
monotonic  increase  of 
response  with  particle 
size  for  both 
nonabsorbing  and 
absorbing  particles. 


Particle  radius  (pm) 

fluctuations  in  the  data  points  where  no  error  bars  are  given,  for  both  non¬ 
absorbing  and  absorbing  particles.  The  theoretical  curves  and  measured  re¬ 
sponse  data  reveal  an  almost  monotonic  increase  in  response  with  increas¬ 
ing  particle  size,  although  the  response  increases  slowly  with  increasing 
size  for  highly  absorbing  particles  with  0.2  /xm  <  r  <  0.6  /xm. 

Collectively,  the  comparison  of  measured  response  with  the  theoretical  re¬ 
sponse  for  both  narrow-angle  and  wide-angle  scattering  probes  reveals 
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that  (1)  the  response  function  (eq  (1))  that  properly  averages  over  intracav¬ 
ity  particle  trajectories  is  in  slightly  better  agreement  with  measurement 
compared  to  the  previous  PAG  response  function  (eq  (3)),  and  (2)  the  fine 
structure  predicted  by  the  response  function  (eq  (1))  is  not  experimentally 
resolved. 

This  latter  finding  is  not  surprising,  since  the  amplitude  of  the  fine  struc¬ 
ture  is  of  the  order  of  the  pulse  height  channel  widths  in  these  probes,  and 
thus  we  should  not  expect  them  to  be  resolved.  Even  the  use  of  a  pulse 
height  analyzer  with  greater  resolution  may  not  resolve  the  fine  structure 
because  of  basic  limitations  of  these  instruments  (laser  intensity  variations 
for  different  particle  trajectories,  etc). 

We  conclude  that  the  theoretical  response  function  (eq  (1))  can  adequately 
predict  the  response  of  both  narrow-angle  and  wide-angle  scattering  probes 
to  spherical  particles,  regardless  of  their  dielectric  properties,  so  long  as 
they  are  of  uniform  composition.  We  provide  such  predictions  for  some 
common  atmospheric  aerosols  in  the  next  section. 


4.  Response  Predictions  for  Homogeneous  Particles 


Atmospheric  aerosols  can  be  irregular  and  of  mixed  composition.  How¬ 
ever,  some  occur  as  spherical  and  homogeneous  particles.  We  consider  four 
types  of  aerosol  particles  that  are  common  in  the  earth's  atmosphere:  wa¬ 
ter  droplets,  droplets  composed  of  a  75-percent  sulfuric  acid /25-percent 
water  mixture,  ammonium  sulfate  particles,  and  particles  of  black  carbon. 
In  addition,  for  generality,  we  consider  metal  particles  of  copper.  We  ini¬ 
tially  assume  all  to  be  spherical  and  homogeneous.  Figures  3  and  4  present 
response  calculations  for  these  aerosols  for  both  narrow-angle  and  wide- 
angle  scattering  probes.  For  reference,  these  calculations  are  compared  to 
the  manufacturer's  calibration,  which  is  plotted  on  the  particle-scattering 
cross-section  scale  after  normalization  by  the  factors  determined  in  the  pre¬ 
vious  section  (from  fitting  the  polystyrene  latex  measurements  to  theory). 


Figure  3.  Predicted 
response  for  PMS 
narrow-angle  scattering 
probes  (ASASP-300  and 
ASASP-300X).  Curves 
for  water,  75%  sulfuric 
acid  in  water, 
ammonium  sulfate, 
black  carbon,  and 
copper  are  shown. 
Response  is  a  sensitive 
function  of  particle 
refractive  index. 
Particles  differing  by  a 
factor  of  two  in  radius 
can  have  same 
response. 
Manufacturer's 
calibration  for  a 
ASASP-300  probe  is 
given  for  reference. 
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Figure  4.  Same  as 
figure  3,  except  for  PMS 
wide-angle  scattering 
probes  (ASASP-X, 
ASASP-100X,  LAS-X, 
LAS-250X,  HS-LAS 
models)  and  also  Royco 
model  236.  As  for 
narrow-angle  scattering 
probes,  particles 
differing  by  a  factor  of 
two  in  radius  can  have 
same  response. 
Manufacturer's 
calibration  is  for  an 
ASASP-X. 


Particle  radius  (|im) 


These  results  confirm  what  is  already  well  known:  the  response  for  both 
types  of  probes  depends  on  particle  refractive  index  as  well  as  size.  Meas¬ 
urement  of  particles  having  a  range  of  compositions  (and  refractive  in¬ 
dexes)  degrades  the  size  resolution  of  the  probes  considerably,  as  particles 
of  radii  differing  by  as  much  as  a  factor  of  two  can  have  the  same  response. 

As  is  evident  from  the  figures,  the  manufacturer's  calibration  passes 
through  the  family  of  curves  that  represent  commonly  occurring  atmos¬ 
pheric  aerosol  constituents.  Thus,  for  measurement  of  unknown  aerosol 
that  consists  of  particles  of  widely  different  composition  (as  for  example  in 
the  earth's  tropospheric  mixed  layer),  the  manufacturer's  calibration,  even 
though  it  divides  size  intervals  into  too  many  channels  to  be  meaningful, 
might  be  appropriate  for  obtaining  average  size  distributions. 

If  the  manufacturer's  calibration  is  used  to  infer  particle  size,  the  errors 
resulting  from  refractive  index  effects  are  more  easily  revealed  in  figures  5 
and  6,  where  we  plot  the  PMS  (manufacturer's  calibration)  indicated  radii 
versus  the  actual  (theoretical  response  function)  radii.  For  both  probes,  the 
PMS-indicated  radius  is  generally  in  error  by  less  than  50  percent.  Note 
that  for  the  narrow-angle  scattering  probes,  the  PMS-indicated  radius  can 
be  multivalued. 
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Particle  radius  (jum) 


These  results  give  a  comprehensive  assessment  of  the  response  of  the  PMS 
probes  to  homogeneous  spherical  particles,  but  what  about  aerosol  parti¬ 
cles  that  are  inhomogeneous? 


5.  Response  Predictions  for  Inhomogeneous  Particles 


For  many  applications,  not  only  does  the  aerosol  of  interest  consist  of  par¬ 
ticles  of  different  composition,  but  single  particles  may  be  composed  of 
mixtures  (black  carbon  in  ammonium  sulfate,  black  carbon  in  quartz,  clay 
minerals  in  quartz,  black  carbon  in  water,  sulfuric  acid  with  crustal  core, 
ammonium  sulfate  in  sulfuric  acid,  silica  shards  in  black  carbon  agglom¬ 
erates,  metals  in  sulfuric  acid,  etc).  Such  particles  can  have  very  compli¬ 
cated  morphology  (Gillette  and  Walker,  1977;  Chylek  et  al,  1981;  Pinnick  et 
al,  1985;  Sheridan  and  Musselman,  1985;  Sheridan,  1989a,  1989b;  Sheridan 
et  al,  1993;  Reitmeijer  and  Janeczek,  1997).  To  predict  the  response  of  the 
probes  to  internal  mixtures,  we  modify  the  response  function  (eq  (1))  by 
choosing  the  simplest  possible  model:  that  of  a  sphere  containing  a  spher¬ 
ical  inclusion  (Ngo  et  al,  1996).  We  consider  four  specific  particle  types:  a 
water  droplet  containing  a  black  carbon  inclusion,  a  sulfuric  acid  droplet 
containing  quartz,  an  ammonium  sulfate  particle  containing  carbon,  and  a 
sulfuric  acid  droplet  containing  metal. 

We  consider  the  narrow-angle  scattering  probes  first.  To  calculate  the  re¬ 
sponse  function  for  black  carbon  in  a  water  droplet  sulfate  host,  we  ne¬ 
glect  the  fractal  cluster  morphology  of  black  carbon  (Forrest  and  Witten, 
1979;  Colbeck  et  al,  1989)  and  assume  a  solid  black  carbon  core  or  inclu¬ 
sion  with  fixed  radius  of  0.1  pm  enclosed  within  the  host.  We  calculate  the 
response  in  two  ways.  First  we  assume  the  core  is  located  in  the  center 
of  the  host  particle  and  increase  the  host  size.  These  results  are  displayed 
by  the  solid  curves  in  figure  7.  Then  we  relax  this  assumption  and  let  the 
inclusion  be  randomly  positioned  within  the  particle,  and  calculate  the  re¬ 
sponse  for  each  position.  The  larger  the  host  particle,  the  more  possible 
positions  are  considered.  The  average  of  these  response  values  is  shown  by 
the  corresponding  open  circles  in  figure  7.  The  standard  deviations  of  the 
response  values  for  different  inclusion  positions  are  also  calculated  but  not 
displayed.  Results  are  calculated  for  only  a  small  number  of  host  particle 
sizes  because  of  the  considerable  computational  requirements.  For  refer¬ 
ence,  the  manufacturer 's  calibration  is  also  shown. 

For  simplicity,  in  calculating  the  response  for  other  internal  mixtures,  we  fix 
the  inclusion  size  to  be  the  same  as  for  black  carbon  in  water  (i.e.,  inclusion 
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Figure  7.  Predictions  of 
response  for  PMS 
narrow-angle  scattering 
probes  (models 
ASASP-300  and 
ASASP-300X)  to 
inhomogeneous 
particles:  black  carbon 
inclusions  in  water, 
quartz  in  sulfuric  acid, 
black  carbon  in 
ammonium  sulfate,  and 
metal  in  sulfuric  acid. 

Solid  curves  are  for  a  ^ 

0.1-pm-radius  inclusion 
at  center  of  particle;  ’•£ 

open  circles  denote  o_ 

average  response  <^L 

values  for  inclusion  o 

positioned  randomly  q 

within  particle.  Data  c 

sets  are  offset  by  1  to  3  q. 

orders  of  magnitude  on  <8 

vertical  scale  as  ^ 

indicated.  Response  of 
a  composite  particle 
rapidly  transitions  from 
that  of  inclusion  (for 
composite  particle 
radius  near  0.1  pm)  to 
that  of  host  as  host 
radius  is  increased. 
Manufacturer's 
calibration  is  shown  by 
solid  dots. 


Host  particle  radius  (pm) 
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radius  r  =  0.1  /im).  The  response  calculations  for  internal  mixtures  using 
this  simple  model  of  a  sphere  containing  a  spherical  inclusion  reveal  the 
following:  (1)  for  small  host  sizes,  where  the  inclusion  nearly  fills  the  host, 
the  response  approximates  that  of  the  inclusion,  but  quickly  transitions,  as 
the  host  size  increases,  to  a  response  approximating  that  of  the  host;  (2)  the 
average  response  for  a  randomly  positioned  inclusion  is  generally  within 
20  percent  of  the  response  for  a  concentric  inclusion;  (3)  the  variation  in 
response  due  to  changing  orientation  is  small,  as  revealed  by  the  standard 
deviations  of  the  response  for  a  randomly  positioned  inclusion  generally 
being  only  about  10  percent  of  the  average  response;  and  (4)  as  for  most 
homogeneous  particles,  the  manufacturer's  calibration  is  not  a  very  good 
approximation  of  the  predicted  response,  since  the  result  in  some  size  re¬ 
gions  is  undersizing,  and  in  others  oversizing. 

Figure  8  presents  the  response  of  the  wide-angle  scattering  probes  to  par¬ 
ticles  containing  inclusions.  Again,  the  average  response  for  a  randomly 
positioned  inclusion  is  generally  within  20  percent  of  the  response  for  a 
concentric  inclusion  (detailed  comparisons  show  that  for  carbon  in  water, 
the  difference  is  less  than  13  percent,  for  silicon  in  sulfuric  acid  3  percent, 
for  carbon  in  ammonium  sulfate  10  percent,  and  for  metal  in  sulfuric  acid 
22  percent).  However,  unlike  for  the  narrow-angle  scattering  probes,  the 
manufacturer's  calibration  is  quite  a  good  representation  of  the  probe's  re¬ 
sponse  for  some  internal  mixtures:  quartz  in  sulfuric  acid,  black  carbon  in 
ammonium  sulfate,  and  metal  in  sulfuric  acid  (so  long  as  the  sulfuric  acid 
mass  fraction  dominates).  The  manufacturer's  calibration  is  less  valid  for  a 
mixture  of  black  carbon  in  water. 
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35°  to  120°  collection 


Figure  8.  Same  as 
figure  7 ,  except  for  PMS 
wide-angle  scattering 
probes  (ASASP-X, 
ASASP-100X,  LAS-X, 
LAS-250X,  HS-LAS 
models)  and  also  Royco 
236. 


Host  particle  radius  (pm) 
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6.  Comparison  of  New  Response  Function  to  PAG  Theory 


What  are  the  differences  between  the  results  of  using  the  new  response 
function  (which  properly  accounts  for  particle  trajectory  through  nodes, 
antinodes,  and  intermediate  regions  of  the  intracavity  laser  beam)  and  the 
results  of  using  the  previously  published  PAG  theory  approximation  (which 
assumes  particle  trajectories  pass  through  only  antinodes  of  the  beam)? 
How  do  these  response  functions  compare  to  the  manufacturer's  calibra¬ 
tion?  To  address  these  questions,  we  have  compared  the  response  functions 
by  making  graphs  of  the  PMS-indicated  size  (manufacturer's  calibration) 
versus  the  actual  size  (predicted). 

Results  for  particles  that  are  nonabsorbing  (m  =  1.50  -  Of)  and  absorbing 
(1.95  -  0.66i)  are  displayed  in  figure  9  (for  narrow-angle  scattering  probes) 
and  figure  10  (wide-angle  scattering  probes).  The  difference  between  the 
new  response  function  (eq  (1))  and  the  PAG  response  function  (eq  (3))  re¬ 
sults  is  small,  particularly  for  highly  absorbing  particles.  For  nonabsorbing 
particles,  the  amplitude  of  the  high-frequency  oscillations  is  much  reduced 
for  the  new  response  function  compared  to  the  PAG  response  function. 
For  highly  absorbing  particles,  the  response  functions  for  the  narrow-angle 
scattering  probes  are  indistinguishable.  A  similar  comparison  for  the  cor¬ 
rected  PAG  response  function  (eq  (2))  in  figure  11  reveals  small  differences 
from  the  original  PAG  response  function  (eq  (3)). 

Figure  12  presents  a  more  direct  comparison  of  the  new  response  function 
and  the  PAG  response  function  for  the  wide-angle  scattering  probes.  For 
nonabsorbing  particles  with  refractive  index  m  -  1.50  -  0/,  the  new  re¬ 
sponse  function  is  a  smoother  function  of  particle  size  compared  to  the  PAG 
response  function. 


19 


Figure  9.  Relation 
between  PMS-indicated 
particle  radii  (smooth 
curve  through 
manufacturer's 
calibration)  and 
predicted  particle  radii 
(new  response  function 
and  PAG  response 
function)  for 
narrow-angle  scattering 
probes.  Solid  straight 
line  is  a  reference 
depicting  perfect 
agreement  between 
indicated  and  predicted 
size.  Curves  for 
nonabsorbing  and 
highly  absorbing 
particles  are  shown. 
Difference  between  new 
response  function  and 
PAG  response  function 
results  is  small  for  both 
absorbing  and 
nonabsorbing  particles. 
Normalization  in  factor 
relating  manufacturer 
calibration  to 
theoretical  response  is 
1  V  =  5.8  X  10-°  cm-. 
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Figure  11.  Relation 
between  PMS-indicated 
particle  radii  (smooth 
curve  through 
manufacturer's 
calibration)  and 
predicted  particle  radii 
(corrected  PAG 
response  function,  eq 
(2),  and  PAG  function, 
eq  (3))  for  wide-angle 
scattering  probes.  Solid 
straight  line  is  a 
reference  depicting 
perfect  agreement 
between  indicated  and 
predicted  size.  Curves 
for  nonabsorbing  and 
highly  absorbing 
particles  are  shown. 
Normalization  factor 
relating  manufacturer 
calibration  to 
theoretical  response  is 
1  V  =  5.3  x  1CT9  cm2. 

Figure  12.  Predictions 
of  response  for  PMS 
wide-angle  scattering 
probes  based  on  new 
response  function  (eq 
(1))  and  PAG  response 
function  (eq  (3)), 
compared  to 
manufacturer's 
calibration.  For 
particles  with  m  =  1.50 
—  0/,  manufacturer's 
calibration  follows  PAG 
theory  and  undersizes 
particles  with  radii  less 
than  0.2  /;m. 
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7.  Predictions  of  Response  for  Multiline  Laser  Operation 


The  manufacturer  claims  (Knollenberg,  1989)  that  the  HeNe  lasers  in  the 
LAS-X  and  HS-LAS  probes  have  sufficient  gain  to  operate  in  multifrequency 
modes,  with  as  many  as  six  lines  from  yellow  to  red  for  the  HS-LAS,  and 
that  the  multifrequency  operation  smooths  the  response  characteristics  of 
the  instrument.  To  address  this  claim,  we  consider  particle  scattering  in  a 
cavity  lasing  in  multiple  lines. 

As  shown  in  appendix  A,  one  can  calculate  the  response  function  for  multi- 
line  operation  by  weighing  the  response  function  for  each  line  by  the  inten¬ 
sity  of  that  line  (eq  (A-2 7)).  According  to  PMS,  the  intensity  distribution  for 
the  HS-LAS  is  typically  80  percent  for  the  632.8-nm  line,  10  percent  for  the 
629.4-nm  line,  and  10  percent  for  all  other  lines  (private  communication, 
John  Mitchell,  Particle  Measuring  Systems,  1998).  We  note  that  by  measur¬ 
ing  the  light  transmitted  through  the  highly  reflecting  closure  mirror  of  a 
factory-refurbished  LAS-X  (serial  number  9075-0786-252,  manufactured  in 
July  1986),  we  found  98-percent  intensity  in  the  632.8-nm  line,  2-percent  in¬ 
tensity  in  the  640.1-nm  line,  and  no  measurable  intensity  in  the  other  lines. 
Thus  the  manufacturer's  claim  of  multiline  operation  may  be  exaggerated. 

Nevertheless,  lasers  in  the  HS-LAS  probes  purportedly  have  higher  gain 
than  those  in  the  LAS-X  probes.  To  assess  the  effect  of  multiline  operation, 
we  assume  the  manufacturer's  claim  is  correct  and  calculate  the  response 
of  the  HS-LAS  (wide-angle  scattering  probe)  for  the  following  distribution 
of  intensities:  80  percent  at  632.8  nm,  10  percent  at  629.4  nm,  and  2  per¬ 
cent  each  at  the  593.9-,  604.6-,  611. 8-,  635.2-,  and  640.1-nm  lines.  Figure  13 
compares  the  results  for  multiline  lasing  and  single-line  lasing  (all  inten¬ 
sities  at  632.8  nm).  Response  calculations  are  presented  for  nonabsorbing 
ammonium  sulfate  and  highly  absorbing  black  carbon.  In  both  cases  the 
difference  in  response  cannot  be  distinguished  in  the  figure  (even  though 
it  may  not  be  clearly  evident,  two  curves  are  displayed  for  each  refractive 
index).  Numerical  results  show  that  the  differences  are  less  than  7  percent 
for  ammonium  sulfate  and  less  than  2  percent  for  black  carbon.  Thus  the 
effect  of  multiline  lasing  on  response  is  negligible. 
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Figure  13.  Predictions 
of  response  for  PMS 
wide-angle  scattering 
probes  having  lasers 
with  sufficiently  high 
gain  for  lasing  in 
multiple  lines  (LAS-X, 
LAS-250X,  HS-LAS 
models).  Response 
curves  are  presented  for 
multiline  and 
single-line  lasing 
operation.  Curves  are 
presented  for 
nonabsorbing 
ammonium  sulfate 
(m  =  1.52  —  Oi)  and 
highly  absorbing  black 
carbon  (m  —  1.95  — 
0.66?)  particles  (two 
curves  are  presented  for 
each  refractive  index). 
Multiline  operation  has 
negligible  effect  on 
response  characteristics. 


Particle  radius  (jjm) 
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8.  Conclusion 


We  present  here  an  improved  response  function  for  light-scattering  aerosol 
counters  that  use  intracavity  laser  illumination.  The  new  response  func¬ 
tion  predicts  a  smoother  size  dependence  compared  to  the  previous  PAG 
response  function,  particularly  for  nonabsorbing  particles,  because  it  ac¬ 
counts  for  phase  averaging.  The  response  function  has  been  used  to  inves¬ 
tigate  the  response  characteristics  of  Particle  Measuring  Systems  aerosol 
probes  that  use  intracavity  laser  illumination.  Probes  having  two  scatter¬ 
ing  geometries  have  been  studied:  those  that  collect  light  scattered  through 
angles  ranging  from  4°  to  22°  from  the  direction  of  the  laser-beam  axis  (as 
in  ASASP-300  and  ASASP-300X  style  probes),  and  those  that  collect  light 
scattered  35°  to  120°  from  the  laser-beam  axis  (as  in  the  ASASP-X,  ASASP- 
100X,  LAS-250X,  LAS-X,  and  HS-LAS  probes).  The  new  response  function 
is  generally  in  good  agreement  with  measurements  on  uniform  particles 
of  polystyrene  latex,  dioctylphthalate,  nigrosin  dye,  and  carbon  black.  Re¬ 
sponse  calculations  for  both  homogeneous  and  inhomogeneous  particles 
with  a  variety  of  refractive  indexes  are  presented  and  compared  to  the 
manufacturer's  calibration;  this  comparison  reveals  that  caution  needs  to 
be  exercised  in  using  the  calibration.  For  the  wide-angle  scattering  probes, 
our  study  suggests  that  the  manufacturer's  calibration  is  generally  appro¬ 
priate  for  sizing  sulfate  particles,  as  well  as  particles  containing  mixtures  of 
black  carbon  in  sulfate  and  quartz  in  sulfuric  acid.  Contrary  to  the  manu¬ 
facturer's  claim,  multiline  laser  operation  has  little  effect  on  probe  response 
characteristics. 


25 


Appendix  A.  Theoretical  Response  Function  for  Intracavity  Laser 
Probes 


A-l  Intracavity  Scattering 

In  active-cavity  scattering  probes,  aerosol  particles  traverse  an  intracavity 
laser  beam  within  an  optical  cell.  Scattered  light  is  collected  either  by  a 
lens  (for  the  narrow-angle  scattering  probes)  or  by  a  parabolic  mirror  (for 
the  wide-angle  scattering  probes).  In  both  instruments,  the  collection  solid 
angle  is  symmetric  about  the  beam  axis,  except  for  small  asymmetry  caused 
by  the  aerosol  entrance  jet  and  exit  port,  which  protrude  into  the  optical  cell 
in  the  X-probes.  These  ports  reduce  in  a  small  but  complicated  fashion  the 
solid  angle  over  which  light  is  collected;  this  effect  is  neglected  here. 

We  consider  the  intracavity  beam  to  be  made  up  of  two  counterpropagat- 
ing  beams.  The  beam  propagating  towards  the  collection  aperture  center  is 
labeled  beam  1  and  the  other  beam  2.  We  define  an  (x,y,  z)  Cartesian  co¬ 
ordinate  system  with  unit  basis  vectors  ( x ,  y,  z )  so  that  the  --axis  extends 
through  the  center  of  the  aperture,  and  z  =  0  is  defined  to  be  the  posi¬ 
tion  where  particles  cross  the  2-axis.  For  our  scattering  calculations,  we  ne¬ 
glect  the  finite  spatial  extent  of  the  beams  and  assume  an  infinite  coherence 
length.  We  model  the  beams  as  linearly  polarized  plane  waves  with  equal 
amplitudes  and  with  parallel  polarization  in  the  .r-direction.  We  use  the 
exp(— itot)  convention  and  represent  the  electrical  fields  of  beam  1  (prop¬ 
agating  in  the  direction  of  increasing  z)  and  beam  2  (propagating  in  the 
direction  of  decreasing  z)  as 

Ej„r(r ,t:  1)  =  xEmcexp[i(kz  -  izt  +  $i)].  (A-la) 

and 

Eindr.t: 2)  =  xEincexp[i{-kz  -  ut  +  $•_>)]•  (A-lb) 

The  wavelength  of  the  beams  in  the  medium  (with  index  of  refraction  N) 
surrounding  the  particle  is  A.  and  the  wavenumber  is  k  =  Nu/c  =  2ir/X. 
We  explicitly  introduce  the  beam  phases  Ty  to  facilitate  discussion  of  laser 
coherence.  The  total  electrical  field  is  the  sum  of  the  fields  of  the  two  beams, 

Einc{r.t)  =  Einc(r,i;p),  (A-2) 

p=1.2 
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and  the  time-averaged  energy  distribution  of  this  monochromatic  field  is 
proportional  to 

Emc(M)  •  E-nc(r,t)  =  |Einc(r,t;  1)|2  +  |Einc(r,f;  2)|2  +  2  Re  [Einc(r,t;  1)  •  E*„c(r,f;  2)] 

=  2(Einc)2[l  +  cos(2  kz  +  3*12)] 

=  4(Einc)2  cos2  (kz  +  $12),  (A-3) 

where  $12  =  $1  —  $2-  We  see  that  maxima  (i.e.,  antinodes  or  bright  fringes) 
occur  in  the  energy  distribution  when  kz  +  $12  =  rnr,  n  =  0,  ±1,  ±2.  ±3, 
and  the  distance  between  the  antinodes  of  this  interference  fringe  pattern 
is<5  =  7r//c  =  A/2.  If  the  coherence  time  of  the  beams  (Tc)  is  longer  than  the 
transit  time  (Tt)  required  for  the  particle  to  pass  through  the  beams,  then 
$12  is  constant  and  the  fringes  are  fixed  in  space,  but  if  Tc  <  Tt,  then  $12 
becomes  a  function  of  time,  as  the  fringes  change  during  transit. 

The  two  incident  electromagnetic  waves  scatter  from  the  particle  crossing 
the  beams,  and  our  immediate  goal  is  to  obtain  the  scattered  wave  solutions 
as  functions  of  the  variables  and  basis  vectors  of  a  single  coordinate  system, 
which  we  choose  to  be  the  (r,  0,  <f>)  spherical  system  (with  unit  basis  vectors 
(r,  0. 4>)  associated  with  the  ( x ,  y,  z)  Cartesian  system.  To  accomplish  this 
goal,  it  is  helpful  to  introduce  two  other  Cartesian  systems:  (x\,y\,z{)  with 
basis  vectors  (x\,yi,z\)  associated  with  beam  1,  and  (:r  >.  y2,  z2)  with  basis 
vectors  ( x2 ,  y2,  £2)  associated  with  beam  2.  The  origin  of  these  two  systems 
is  the  center  of  the  particle,  and 


x2  =  x , 

-m  =  y,  (A-4) 

-z2  =  z. 

As  the  particle  crosses  the  beams,  the  large  spatial  dimension  of  the  col¬ 
lection  aperture  compared  to  that  of  the  aerosol  stream  that  crosses  the 
beam  (about  100  /tm)  and  particle  diameter  (a  few  micrometers)  allows  us 
to  make  the  approximations 


x\  = 

y\  = 

£1  = 


X\  =  x2  =  X. 

yi  =  -2/2  =  y-  (A-5) 

Z\  —  — 2  2  —  Z. 


28 


The  spherical  systems  associated  with  these  two  Cartesian  systems  are 


=  (r,9,d>), 

(r.  0-2-,  62)  =  (r,  7 x-9,~4>),  (A-6) 

with  unit  basis  vectors 

(r,9  i,0i)  =  ( r,9,4> ). 

{r,9-2,4>2)  =  (r,-9.-4>)-  (A-7) 

We  now  write  the  electrical  fields  of  the  two  incident  waves  (p  =  1,2)  as 

Einc(r ,t:p)  =  xpEinc  exp [i{kzp  -ut  +  $p)],  (A-8) 

and  the  far-field  scattering  solutions  as 

ESCfl(r  .t\p)  =  Esca(r;p)  exp  [i(-upt  +  $p)],  (A-9a) 

where 

Esca  (r;  p)  =  (E0p0p  +  E^p)  .  (A-9b) 


The  scattered  fields  can  be  expressed  as  functions  of  amplitude  scattering 
matrices  (Bohren  and  Huffman,  1983)  with  components  Sk  (k  =  1  to  4),  and 
for  incident  fields  polarized  in  the  xp  direction. 


’  Eep  ‘ 

E  ( exp(ifcpr)\ 

Sziftp) 

1 - 

3 

^0 

COS  (f)p 

,  (A-9c) 

m  . 

inc\  —ikpr  J 

_  5*4  (/ip) 

Si(fJ>p)  J 

1 - 

£ 

1 _ 

where  pp  =  cos  9P  and  kp  =  Nup/c.  If  the  particle  is  a  homogeneous  sphere, 
then  S3  =  S4  =  0. 

Recalling  equations  (A-6)  and  (A-7),  we  find  that 

ESCa(r;  1)  =  Edi9  +  E^4>. 

Esco(r;  2)  =  -Ee.29  -  E^l.  (A-10) 

If  the  particle  does  not  pass  perpendicular  to  the  beam,  then  it  has  a  small 
velocity  component  parallel  to  the  beam  axis,  and  the  frequencies  of  the  two 
scattered  waves  are  slightly  unequal  because  of  small  Doppler  shifts.  An 
appropriate  expression  given  by  Jackson  (1975,  eq  (11.8))  can  be  rewritten 
to  obtain  the  frequency  of  incident  beam  p  in  the  particle's  rest  frame  (PRF) 
as 

(ujp)prf  =  oj  —  kv  ■  zp.  (A-lla) 
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and  a  second  Doppler  shift  of  scattered  wave  p  into  the  lab  system  gives 
the  frequency  of  the  far-field  scattered  wave  at  a  point  on  the  aperture  (in 
direction  f)  as 


UJp  —  ( UJp)pRF  —  fcv  •  f, 

—  oj  -J-  Acj^,  (A-llb) 

where 

Ac Op  =  —  kv  •  (zp  +  f),  (A-llc) 

so  that  the  difference  between  the  scattered  frequencies  at  a  point  on  the 
aperture  is 


UJ  21  =  C^2  —  C^i , 

=  2fc(v  •  £), 

=  2fcu2.  (A-12) 


The  Doppler  shift  of  each  beam  is  position  dependent  over  the  aperture,  but 
the  difference  in  frequencies  between  the  two  scattered  waves  (differential 
Doppler  shift)  is  not  position  dependent. 

The  total  electric  and  magnetic  fields  at  point  r  on  the  aperture  at  time  t  are 

Esca(r,f)  =  y  Esca(r ,t\p). 
p= 1,2 

Hsca(ri)  =  Y,  Hsca(r,f:p).  (A-13) 

p=  1,2 

and  the  instantaneous  Poynting  vector  is 

S(r .#)  =  Re[Esca (r i)]  x  Re[Hsca  (r.f)].  (A-14) 


We  assume  that  all  radiant  power  scattered  into  the  aperture  solid  angle, 

PAn(0=  f  dnr2r-S(r.t).  (A-15) 

J  aq 

is  delivered  to  the  sensing  element,  which  produces  an  electronic  signal 
proportional  to  this  power.  We  find  that 


f-s(r’*)  =  ^ReE  E je. 

r  p=l,2  q=l,2 


(tip)  ■  EsC0(r;  q)}  exp [i(-upqt  +  $pq)}, 

(A-16) 
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where  upq  =  up  -  tcq  and  $pq  =  The  time-averaged  irradiance  of 

each  incident  beam  is 

line  =  7T-( Einc f,  (A-17) 

2fj,c 

so  that 

PAtt(t)/Iinc  =  Re  ^  ^  Mpg  exp^u^f  +  %q)],  (A-18a) 

p=  1,2  9=1,2 


where 

Wp<i  =TT  [  dn{kpr){kqr)Esca(r;p)-El^(r;q)/{Einc)‘2.  (A-18b) 
Kpkq  J  Af2 

We  observe  that  W*q  =  Wqp  so  that  W\\  and  11' 22  are  real. 

In  order  to  obtain  a  realistic  simulation  of  the  response  function,  we  need 
to  consider  the  averaging  effect  of  signal  processing  by  the  probe  sensor 
electronics.  To  do  this  we  examine  the  passage  of  the  particle  through  the 
beams  from  entrance  at  time  t  =  0  until  exit  at  t  =  Tt.  After  the  particle 
enters  the  beams,  scattered  light  begins  to  reach  the  sensor,  producing  a 
signal  that  peaks  as  the  particle  reaches  the  maximum  irradiance  at  the  z- 
axis  when  t  ~  Ttj 2.  The  instrument  response  to  a  particle  passing  through 
the  beams  is  obtained  by  electronic  averaging  of  the  signal  generated  by 
the  sensor  over  an  interval  of  time  Te  (a  characteristic  RC  time  constant 
of  the  sensor  amplifiers)  centered  on  the  occurrence  of  the  peak  signal  at 
t  ~  T)/2.  We  compute  the  instrument  response  R  by  averaging  equation 
(A-18a)  over  the  appropriate  interval  of  time. 


R  —  (T>Afl(f  )/Iinc) 

l  rt=T,/2+Tt/2 

=  —  dt[PAn(t)/Iinc], 

Te  Jt=Tt/2—Te/2 

=  Re  5Z  5Z  +  $pq)}) 

p—1.2  q— 1.2 

=  Ik'n  +  W22  +  2R.e[ITTi2(exp[i(u;2i#  +  $12)])]- 


(A-19) 


Thus  the  response  function  consists  of  three  cross-section  components:  one 
associated  with  beam  1,  one  with  beam  2,  and  an  interference  term. 
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A-2  Corrected  PAG  Response  Function 

We  now  derive  a  new  expression  of  the  PAG  theory,  which  is  applicable 
when  the  particle's  trajectory  is  exactly  transverse  to  the  beams  (i.e.,  u>2i  =  0 
and  ki  =  k2  =  k)  at  the  position  of  an  unmoving  bright  fringe  (antinode) 
so  that  $12  =  0.  These  assumptions  reduce  equation  (A-19)  to 

R  =  Wu  +  W22  +  W12  +  W2i  ,  (A-20) 


where 


Wn 
TT22 
Wi2 
IV 21 


j[  min 

¥ 


[ . dp  ||Si(/,)|2  +  |52(m)|2  +  |53(m)|2  +  |S4(/i)| 

J  (A max 
7P  min 


£  r  ^  [l5i(-^)l2  +  ^(-^l2  +  |S3(-/r)|2  +  |5'4(-m)|2]  , 

^  J  ft  max 

-2  /  d/i  [S'1(/i)5*(— /x)  -  52(/i)55(-/i)  +  S3(m)^(-/i)  -  54(/i)54*(-/i)]  , 

^max 
7T  /'Mmin 

T2  /  dii  [S^Sii-ii)  -  SZ(ii)S2(-/i)  +  SI(ii)S3(-fi)  -  Sl(/i)S4(-ii)} ; 

^  J  fJ> max 


and  where  /imin  =  cos  #mjn  and  //max  =  cos  #max .  We  observe  that  equation 
(A-20)  can  be  rewritten  as 

R  *  4  r*  [|5i(/i)  +  51(-/i)|2-H52(/x)-52(-M)|2 

^  ^Mmax  L 

+  |53(/i)  +  53(-M)  I2  +  |54(/i)  -  S4(-/*)|2]  ,  (A-21) 

and  if  the  particle  is  a  homogeneous  sphere,  then  the  integrands  S3(/t)  = 
54(/t)  =  S3(— /t)  =  S4(— /1)  =  0.  and  we  obtain 

f?  -  TrAr2  rmin  d/x  [  |5i(/x)  +  Si(-/x)|2  +  |52(/x)  -  52(-/x)|2]  .  (A-22) 

J  M  max 

We  refer  to  this  result  as  the  corrected  PAG  equation,  because  it  is  a  cor¬ 
rected  version  of  the  original  PAG  equation  given  by  equation  (3)  in  the 
main  body  of  the  report. 


A-3  New  Response  Function 

To  obtain  a  more  realistic  response  function,  we  now  consider  particles 
crossing  the  nodes,  antinodes,  and  intermediate  regions  during  transit 
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through  the  intracavity  beam.  If  the  beam  fringes  do  not  change  or  move 
during  particle  transit,  then  $12  remains  constant,  and  the  interference  term 
in  equation  (A-19)  can  be  simplified,  giving 

R  =  Wn  +  W22  +  2{sin(w2iTe/2)/(w2iTe/2)}  R*{Wi2  exp[i(u2iTt/2  +  $12)]}-  (A-23) 

Recalling  equation  (A-12),  we  find  that  uj21T€/2  =  2ir(veTe/X).  To  estimate 
the  magnitude  of  this  term,  we  assume  the  diameter  of  the  aerosol  stream 
that  intersects  the  laser  beam  to  be  100  /  /  m ,  and  the  particle's  velocity  trans¬ 
verse  to  the  beam  to  be  10  m/s,  so  that  the  transit  time  across  the  beam,  T%,  is 
10  [is.  Also,  we  assume  that  the  particle's  velocity  is  not  exactly  transverse 
to  the  2-axis  but  has  a  component  v~  that  is,  let  us  say,  2°  from  the  perpen¬ 
dicular  direction,  so  that  (u2)max  =  (10  m/ s)  tan  (2°)  =  0.349  [im/ [is.  Using 
A  =  0.6328  yum  and  Te  =  1.0  [is,  we  find  u2iTe/2  =  2tt  (0.349)(1.0)/0.6328 
and  sin(w2iTe/2)/(w2iTe/2)  =  -0.092.  Thus  this  factor  causes  a  91-percent 
reduction  of  the  magnitude  of  the  third  (interference)  term  in  equation  (A- 
23).  This  reduction  is  even  more  for  larger  values  of  ui2iTe/2  =  n(vzTe/5). 
This  last  expression  ofu>2iTe/2  is  equal  to  tt  times  the  number  of  fringes  that 
the  particle  crosses  during  the  time  Te,  so  the  interference  term  of  equation 
(A-23)  becomes  smaller  if  the  particle  crosses  more  fringes. 

Having  demonstrated  that  the  interference  term  of  equation  (A-23)  is  small 
for  realistic  particle  trajectories  through  the  beam,  we  turn  our  attention  to 
an  additional  effect  that  may  also  reduce  the  magnitude  of  this  term:  that 
of  finite  laser  coherence  time.  We  examine  the  hypothetical  situation  when 
the  particle  passes  exactly  perpendicular  through  the  beam  so  that  vz  =  0 
and  u>2i  =  0  (i.e.,  there  is  no  differential  Doppler  shift).  Even  though  there 
is  no  Doppler  shift,  the  fringes  are  not  fixed  during  particle  transit,  because 
the  laser  changes  longitudinal  modes  ($12  is  not  constant,  as  the  coherence 
time  of  the  laser,  Tc  =  0.5  [is,  is  less  than  the  particle  transit  time  and  the 
electronic  averaging  time).  Thus,  $1  and  $2  each  assume  multiple  random 
values  during  particle  transit,  and  we  find  that 

1  ft=Tt/2-\-Tefo 

(exp[i(u;2i*  +  $12)])  =  7f~  /  dtexp[i$i2{t)]  —  0.  (A-24) 

J-e  Jt=Tt/2-Te/2 

We  conclude  that  the  response  function  that  accounts  for  particles  crossing 
fringes  of  the  intracavity  beam  and  the  finite  coherence  time  of  the  laser  can 
be  approximated  by 
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R  =  Wn  +  W22 

77-  fV>  min  r  r»  o  o  o 

=  Z2  /  ^  |5i (p)|2  +  l^2(/u)|2  +  |53(m)|2  +  |54(m)|2 

^  ''/Umax 

+  |51(-/i)|2  +  |5'2(-p)|2  +  |S3(~p)|2  +  |S4(-p)|2] ,  (A-25) 

where  we  have  used  the  simplifying  approximation  that,  since  particle  tra¬ 
jectory  velocities  are  within  a  few  degrees  of  being  perpendicular  to  the 
beam,  l/(k\)2  cs  1  /(A^)2  —  l/fcifo  —  1  /k2. 

If  the  particle  is  a  homogeneous  sphere,  then  the  integrands  S3 (p)  =  S4(p)  = 
S3(— p)  =  S4(— n)  —  0,  and  we  obtain  the  new  response  function: 

r  *  Y2  r*  d^i(^)i2 + i5'2(>u)i2 

™  ^Mmax 

+  l'S'i(— p)|2  +  |52(— p)|2J  .  (A-26) 

A-4  Response  Function  for  Multiline  Lasing 

Finally,  we  consider  the  instrument  response  for  multiline  lasing.  If  the  in¬ 
tracavity  beam  contains  multiple  frequencies,  then  the  difference  between 
these  frequencies  is  approximately  a  factor  of  106  greater  than  the  differen¬ 
tial  Doppler  frequency  considered  earlier,  and  the  time-dependent  interfer¬ 
ence  terms  in  the  total  instantaneous  Poynting  vector  are  averaged  to  zero 
by  the  sensor.  The  result  is  that  the  time-averaged  instantaneous  Poynting 
vector  is  equal  to  the  sum  of  the  time-averaged  Poynting  vectors  of  the 
different  frequencies,  and  the  instrument's  response  is 

E(pAn(i,A)) 

R  =  -*= - 

/  /inc(^) 

A 

J>(A)J<nc(A) 

_A _ 

Ewa) 

A 

=  £H(A)/(A),  (A-27a) 

A 

where  /(A)  is  the  fraction  of  the  total  irradiance, 

/(A)  =  Iinc(X)/ E  W A).  (A-27b) 

A 
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Thus  the  multifrequency  response  is  the  weighted  sum  of  the  response 
functions  (of  wavelength  A),  each  weighted  by  the  irradiance  fraction  at 
A. 
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Appendix  B.  Normalization  of  Measured  Response  to  Theory 


The  measured  responses  of  the  PMS  probes  are  fit  to  the  theoretical  re¬ 
sponse  as  follows.  Measurements  of  the  raw  response  of  a  particular  probe 
(in  volts),  Rraw(rn):  where  r„  is  the  polystyrene  latex  particle  radius,  are 
taken  from  the  appropriate  journal  reference  for  a  total  of  N  particle  sizes. 
We  restrict  the  fit  to  the  polystyrene  latex  data  because  these  particles  have 
well-characterized  size  and  refractive  index.  The  experimental  response  of 
the  instrument  (in  centimeters  squared),  R^p,  is  then  given  by 

R™p  =  GRrau’(rn),  (B-l) 

where  G  is  the  gain  (in  centimeters  squared  per  volt)  to  be  determined  by 
comparison  of  the  experimental  response  with  a  theoretical  cross  section 
(in  centimeters  squared).  The  corresponding  theoretical  response  Rsymffn) 
is  computed  with  equations  (1),  (2),  or  (3)  in  the  main  body  of  the  report. 

The  approach  used  here  is  to  determine  the  gain  G  by  minimizing  the  ab¬ 
solute  value  of  a  function  S  (G)  defined  as 

S(G)  =  ^2  ( [kJg  10 ( R^n  ^ ) ]  -  P°§  1 0 (R-sym ( ) ] )  •  (B"2) 

n=  1 


This  approach,  in  which  the  sum  of  the  difference  of  the  logarithm  of  re¬ 
sponse  is  minimized,  weighs  all  response  measurements  nearly  equally  in 
the  fitting  procedure. 
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